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Tbermocouples (t/c) positioned in the furnace element and 
process chamber of an oxidation/diffusion or LPCVD 
furnace system are used as inputs for adjusting the power 
d livered to muftizone heating elements. Until recently, the 
industry relied on proportional integral derivative (PID) 
algorithms to indirectly control the temperative of the 
wafer based upon input from the thermocouples. 

Model-based temperature control (MBTC) uses a math- 
ematical model of the thermal characteristics of the fur- 
nace, process chamber, and wafer stack to predict and 
control the wafer's temperature. This article descrit>es the 
implementation, testing, and results of MBTC on a vertical 
thermal reactor system. Significant t>enefits result from 
improved uniformity, rornto-run repeatability, reliability, 
and cycle-time reduction for selected oxidation/diffusion 
processes. 

In the fabrication of semiconductor de\nces, oxidation and dif- 
hision furnaces elevate the silicon substrates to a sufficient tem- 
perature to carry out diffusion of a dopant or oxidation on the 
surface. T\'pically, oxidation and diffusion furnaces are config- 
ured in a similar manner with the same reactant gas and tem- 
perature capabilities, allov^ing most of these sj'Stems to perform 
the same tasks in a production line. Process temperatures range 
from4avi25Cy='C 

For a desired oxidation or diffusion process, wafers are loaded 
into a quartz fixture called a "boat/' as shown in Fig. 1. The 
boat is then inserted or "pushed" slowly into the furnace 
tube which is controlled to a temperature above room tem- 
perature but low enough to prevent thermal stress damage 
to the wafers during the insertion. Push temperatures in this 
paper range from 650-800''C. After boat insertion, the furnace 
temperature is ramped to the desired diffusion or oxidation 
temperature, and allowed to stabilize; then the oxidation or 
diffusion step is executed in an appropriate gas ambient. An 
inert gas purge or anneal sequence then occurs, typically at the 
same temperature, followed by a temp>erature ramp down to 
the desired wafer extraction or "pull" temperature. The pull is 
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ngur* 1. Schematic of a typical vertical thermal reactor (VTR) showing 
the location ot spike and profile thermocouples. Neither set directly mea- 
sures the wafers, and the discrepancy t}€tween measured ar)d wafer tem- 
peratures results in reduced process control. 

normally performed at the same temperature as the push, and 
the wafers are slowly removed from the furnace tube into 
the room temperature environment. 

The general reactor configuration includes a process chamber 
made up of a high-purity quartz process tube, wafer boat, and 
pedestal, which the wafer boat resides on (Fig. 1). A hydrogen 
combustion system is housed in an externally configured heating 
element. The process chamber is positioned vertically inside a he- 
lically wound healer element with five control zones. A SiC 
liner is positioned between the process chamber and the heater 
element. Gas injection is accomplished \ia mass flow controllers 
(MFCs) through a manifold assembly that directs process gases 
into the top of the process tube. 

As de\ice geometries continue to shrink, chip manufacturers 
must achieve proportional improvements in the control of gate 
oxide unifonnities, line widths, diffused junction depths, and 
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related dimensions. Comparable technology advances must be 
realized in difhision furnace temperature control. 

Furnace temperature is monitored and controlled via ther- 
mocouples located strategically to represent each of the five control 
zonesof the heater element "Spike" thenmooouples are positioned 
outside of the process chainba' in dose proximity to the dement 
windings in each of the five c nlrol zones. Readings from these 
spike thermocouples approximate the heater elonent tempera- 
ture. "ProfiJe'' thermocouples are positioned inside the process 
duunber in dose proximity to the wafer stack, one in eadi of ^ 
five control zones. These profile thermocouples approximate 
the wafier stack temperature. 

A vertical diffusion funuice may process as many as 150 pro- 
duction wafers in a single batch. To meet the needs of next geiv 
eration device geometries^ the entire batch must have tighter 
thermal uniformity and lower expenditure of thermal budget 
during ramps and stabilization. KfBTC technology successfully 
addresses these issues and provides the wafer temperature control 
required to meet present and future batch manijfacturing needs 
of the semiconductor industry. This artide describes the results 
of a joint development agreement QDA) between Silicon VaUey 
Group's (SVG) Thermco Division and Motorola's COM-1 to im- 
plement MBTC on a vertical fumace. 

Evaluation metliocfs 

Two metrics were used to evaluate the performance of temper- 
ature control systems: 

1 ) direct measurement of wafer temperature or "thermal uni- 
formity" and 

2) oxide thickness uniformity. 
Thermal uniformity results were obtained using type R ther- 
mocouple instrumented wafers from SensArray. Each wafer had 
a thmnocoupie located at the wafer center and the wafer edge. In 
addition, profile and spike temperature data was collected on 
all runs. Wafer thermocouple data was fed into a validated two- 
stream oxide growth model to compute variations in oxide thick- 
ness for each wafer thermocouple location. This model was 
used to compute the target oxide thickness. Da, using the setpoint 
temperature, arni actual oxide thickness, Di, using measured wafer 
temperatures. Using these quantities, two measures of oxide uiti- 
fomiity were computed: 
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with respect to the value for the temperature set-point trajec- 
tory. Thus, 
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where 

n = the number of measurement locations considered 
s the mean of all D,. 

Computed oxide variations were converted to equivalent tem- 
perature variations using sensitivity of oxide thickness to tem- 
perature, Sttt (A/°C), computed from the oxide growth model. 
Thermal uniformity is computed two ways. Self-imiformity is 
computed with respect to its own mean. Total thermal uniformity 
(TTU) is computed for all wafer thermocouple locations in a batch 



Self thermal uniformity = (^C) (2b) 

TTU values are larger than self-uniformity numbers. Smaller 
TTU values ensure that run-to-run variations and fumace-to- 
fumace variations are smaD. 

• Oxide titickness imiformity tests were f>erfmned using lightly 
boron doped <100> crystal orientation test wafers. All pro- 
cessing was performed using 150-mm wafers and evaluated 
on an ellipsometer using a nine-point measurement pattern with 
3-mm edge exclusion. Oxide uniformity results are reported 
as % iKi values calculated by dividing 1-a A x'aluesby theaverage 
thickness x^. 
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The wafer load for all funs was based on the Motorola 
COM-1 production load size of 37 wafers with a wafer spadng of 
0317 in. The thermal uniformity runs had 5 t/c wafers placed 
in the load at wafer positions 6, 12, 19, 26, and 32 (counting from 
the top down); the rest of the load was made up of filler waf«3. 
The t/c wafer runs were performed in a nitrogen ambient lo 
prevent oxidation of the SensAiray wafers. The t» uniformity runs 
had test wafers loaded in positions 6, 19, and 32. 

Proeess descriptions 

Three Motorola oxidation and diffusion p>rocesses were used to 
evaluate and compare PID \'s. MBTC. These processes were chosen 
due to their inherent sei^sitivity to temperature control, which 
greatly influences oxide thickness and diffused junction depth 
uniformity. The process ledpes woe Iransferred to the lab funwce 
and implemented in both PID and MBTC temperature control 
versions. The 200 A and 350 A processes were standard oxidati n 
cydes whidi ramp up to oxidation temperature and stabilize in 
an inert gas, oxidize in either dry O2 or steam, perform a post- 
oxidation aiuieal, and ramp down to the pull temperature in an 
inert gas. All processing was performed at atmospheric pressure. 
A production emitter anneal process is normally performed en- 
tirely in an inert gas to diffuse the emitter junction to a desired and 
repeatable dq>th within the base region. This process was mod- 
ified for this project by replacing the inert gas uith Q2 throughout 
the entire process, including the push and pull. In this way, the 
time-averaged thermal uniformities could be evaluated by mea- 
suring oxide thidcnesses. 

PID eontrol limitations 

The industry standard PID control algorithm adds three terms to- 
gether to determine power applied to each zone of the heater 
element. Each term is a function of only the error between the 
measuremait and the setpoint for the control loop. By design, a 
PID controUo- can only act on tempmture measurements local 

contmuea on page 122 



1 20 Sobd State Technology Nemmber 1 998 



Furnace SdBTC continued from page 120 



TaM 1p Comparison of PID ai 
MBTC parf Ofmanca results 
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to the zone being controlled by it. 



Error e = Setpoint - Measurement 
Heater power = K*e t- K^ ^de/dt + Kj* jedt (4) 

The (first) proportional term acts as a primary negative feed- 
back on the setpoint error. The (second) derivative term acts as 
a brake to reduce speed of response, and the (third) integral 
term acts to reduce steady-state offset in the error. 

PID contrdlers have t>een popular o\'er the last 40 years because 
they could be implemented using analog components. When mi- 
croprocessors came along, digital PIDs were simple enough to be * 
coded and implemented viith the modest computation resources 
available on S-bit processors. PIDs are tuned using riles of thumb, 
such as the Zielgler-Nichols rules, directly from observed thermal 
step responses on the furnace. 

Spike, Profile, Ratio, and Cascade control strategies are used on 
production furnaces. Spike and Profile strategies use the spike and 
profile t/c measurements in forming the setpoint error. Ratio 
control uses spike and profile t/c signals in a fixed ratio to form 
the measurement signal Cascade control refezs to a sehes or cascade 
of two PID loops per zone. In addition, programmable limits 
are placed on power outputs. 

Tuning rules for PID assume there is no coupling among heating 
zones. However, the profile thermocouple in the end-zone is af- 
fected by power applied to the adjacent zone. This coupling makes 
tuning difficult and places inherent limitations on performance. 
For Qcample, tightening cascade ODntrol or profile control to reduce 
offsets leads to oscillations. A spike thermocouple sees less 
effect of power in adjacent zones than the profile thermocouple 
because the spike thermocouple is closer to the heater element. 
Therefore, spUce PID control can be tuned to respond much faster 
than profile or cascade control. Spike control is used during sta- 
bilization and processing to ensure fast PID response. 

After profiling, the steady-state offset between profile and 
spike readings is estimated and stored in a profile table. This 
offset is used to bring steady-state profile tempjeratures to the 
processing temperature by adding a feed-forward bias to the 
spike control loop. 

In practice, furnace reprofiling takes up to 12 hours and can 
be required as often as oiwre a week. In some processes, where 
cycle time and thermal budget pennit, a real-time profiling step 
of up to 30 min is added to the recipe at the process temperature. 
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PID performance baseline 

Process perf ormanoe relative to the current Motorola performance 
results was verified using standard PID temperature contrcsl. Eadx 
PID process was then baselincd by performing passive data col- 
lections (PDCs) consisting of five consecutive runs each. Passive 
data collection is performed \\ith a strict "hands off" approach 
where no ad justnients are made to the process or system through- 
out aU runs. The tox uniformity and throughput performance results 
for the PID baseline are depicted in Table 1 . The % 1 -o is calcu- 
lated using all points from the five runs, with each PDC having 
135 data points. 

figure 2 plots temperature deltas for the Sens Array t/c wafers 
vs. the profile thermocouples. This plot displays temperature delta 
vs. time from the start of ramp up through the oxidation or anneal 
step to the end of the ramp-down step. Two main observations 
are evident: 

1) profile t/c temperature is not an accurate approximation of 
wafer temperatures, and 

2) PlD temperature control does not take into account wafer 
temperature. Even with finely tuned PID, wafer to profile t/c 
temperature deltas will always be present, since true wafer 
temperatures are unknown to PID control. 

The largest deltas, up to 15**C, exist during the ramp up (from 
42 to 62 min) and ramp down (from 95 min to the end of the 
chart). During the oxidation or anneal (steady-state) step the 
wafer tem]:>eratures are within 5**C of the respective profile ther- 
mocouples (Fig. 3). For this process, PID shows good control 

continued on page 124 




Time (miii) 



Figure 2. The difference between pfofile thermocouple and wafer edge 
thernDocoupie lor each zone shows large differences during the ramp-up 
and ramp-down (dynamic) measurements. 




Figure 3. The a) profile and b) wafer-edge temperatures show good control 
under PID for the anneal process. The profile temperature undershoot at tt>e 
end of the ramp up keeps the wafer temperature from overshooting. 
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Furnace MBTC continued from page 122 

of the wafer lemperntures to the temperature setpoint. The under- 
shoot at the end of the ramp keeps the wafer temperature 
from ovcrshofUinf;. 

MBTC furnace configuration 

MBTC ^\'as implemented on a vertical thermal reactor (VTR) model 
\nrR7000 l.ocated in the Thermco Division's applications bbora- 
tor>' in San Jose, CA, as well as on a production VTR70CK) located 
in Motorola's COM-1 facility in Phoenix, AZ. All de\'elopment 
work and perfomiance e\'ahiation was perfbrmed on the SVG ap 
plications lab VTR. The process chamber and element configu- 
ration of the XTR in the lab is a duplicate of the VTR configura* 
tion at COM-1, facilitating a smooth MBTC technology transfer 
to the COM-1 \TR. 

The term Model-Based Control is used to denote a broad 
class of control approaches using systematic process modeling to 
determine control algorithms. In the past, practical applications 
t)f M(Kiel Based Control have been limited to control of equipment 
\'ariables measured on-line, such as the profile and spike ther- 
m<xrc^upJes. 

The MBTC approach described here is unique in the use of 
mixiels to create production-worthy virtual smsom for wafer tem- 
peratures and silicon parameters and in the use of these sensors 
for real-time control. RdMan's Thermal Reactor Optimizer-1 (TRO- 
1) software uses an advanced mathematical thermal model of 
wafer^i and the wafer en\Tronment (furnace/process chamber) to 
estimate and control within wafer and wafer-to-wafer tempera- 
tures. It uses existing spike and profile thermocouples as well as 
heater powers to control heaters in real-time. 

In order to implement MBTC on the VTR, several hardware 
«ind software changes were required. These included upgrad- 



ing the CPU used for temperature control, inoorpomting soft%vaTe 
to perform MBTC. adding hardware to measure heating ctemenl 
power in each zone, and upgrading the user-into-fece soft^^are. 
One of the requirements of the MBTC algorithm is measure- 
ment of the actual power being delivered to the element. The stan- 
dard VTR configuration does not have funvKe elem nt power 
measurement capabilities. Additional hardware was thus desigrted 
into the VTR to sci^ current and voltage for each zone and report 
a power value back to the MBTC CPU board. MBTC param- 
eters were added to the data collection capabilities and furnace 
configuration menus. 

Oxidation control 

In most oxidation processes, wafer tem p er a ture adherence to 
the process setpoint is most important during the oxidation 
step. In particular, for optimum oxide thickness uniformities and 
oxide integrity, it is more critical to control wafer temperature 
during and iirunediately after the oxidatk>n than during the ramp 
up and ramp down. In this type of process we incorporate a MBTC 
strategy called flat trajectory, in which the model drives the wafer 
temperatures to the process setpoint and hoWs them at setpoint 
for a specified time. The objective is to tightly control wafer 
temperatures during oxidation and armeal for the full batch of 
processed wafers. In submicron oxidation processes, reducing 
ramp and stabilization time significantly reduces expenditure 
of thermal budgets. 

Diffusion control 

For a diffuson process, wafer temperatures must be precisely con- 
troUed through the entire thermal cyde. Total thermal uniformity 

contatuea on page T26 
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Furnace MBTC continued from page 124 




n^ure 4. Calculated oxide growth as a function of time for each wafer 
thermocouple location during the anneal process: a) MBTC oxide growth 
trajectory, and b) PID oxide growth trajectory. The MBTC dynamic optimiza- 
tion strategy minimizes variation at the end of the process, resulting in an * 
order of nfiagnitude improvement in TTU over baseline PID. 



(TTU) of the entire batch during ramp up, stabilization, and ramp 
down is critical to the control of junction depths and channel 
widths. Using a dynamic optimization MBTC strategy, tlie model 
drives the wafers to the appropriate temperatures to reach a 
low TTU value such that, at the end of the dif^ion cycle, all wafers 
have been subjected to the same thermal budget. The model adjusts 
and optimizes for temperature nonuniformities across each wafer 
and from wafer-to- wafer in the boat. For example, the bottom* ' 
most wafers in the stack may be nominally cooler than the top- 
most wafers during the ramp up. The model achieves good 




1 



Fltfur* 5. Wafer temperatures for a) model-based 200-A oxidation and 
b) PID 2(X)-A oxidation. 
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TTU by drixing the bottom wafers sufficiently hotter than the top 
wafers at the end of the ramp cycle and into'the drive-in step. 

MBTC process optimization and results 

Both modeling approaches use furnace fingeqpiinting to pro\ide 
the empirical inputs used in achie\nng optimal TTU. The finger- 
printing operation serves two important objectives. First, it 
characterizes the unique thermal behaxnor of the furnace config- 
uration. Second, it verifies the integrity of all thermal subsys- 
tems of the furnace. For fingerprinting, wafer mounted thermo- 
couples, such as SensArray wafers, are processed through the 
furnace in an inert atmosphere. To provide a precise and accurate 
model die finger printing data must indude heater powers, wafer 
temperatures, profile temperatures, spike temperatures, and 
gas flows. 




ngure 6. Wafer temperatures for a) model-based DO emitter anneal 
process, and b) PIO emitter anneal process. 



A specially designed fingerprinting kit is used during final test 
at the factory which includes the quartz boat and SensAnay wafers. 
Using the fingerprinting data, a customized mathematical algo- 
rithm is created for each furnace. This allows process engineers to 
establish recipes that are relatively independent of different tcx>ls. 
Because of unavoidable differences among fumaces, a prcxJuc- 
tion PID amtrollcr requires furnace-specific adjustments to recipe 
setpoints. By eliminating these adjustments, MBTC makes rodpe 
transfer easier. MBTC also eliminates the profiling that is required 
with PID control. 

The MBTC x'ersion of each prcKiess was optimized ti> achie\'e the 
desired goals. Tlie flat trajecton' model was used for the two oxi- 
dation processes because the gtvjl was to maximize thn>ughput. 
Temperature ramp rates could be increased and stabilizations 
reduced or eliminated. R>r the 2C>0 A dr>- oxidaticm, the ramp-up 
and stabilization time was redviced from 68 to 50 min. Tlie ramp- 
down time was also reduced from 11 7 to 1 02 min for a total 
cycle-time reduction of 33 nun. For the 350-A wet oxidation, the 
ramp- up and stabilization time was reduced from 44 to 25 min. 

conttnuea on page 12S 
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Furnace M07C continued from page 126 




Fitfnro 7. Water temperatures for a) model-based 350-A oxidation, and 
b) PID3S0-A oxidation. 




1 2 3 4 5 6 7 8 : 9 10 It 




335 I , 

123456 789 10 11 
Run ID 



Figtfra 8. Comparison of predicted oxide mean thicknesses calculated 
from water temperature measurements made during the Phase I Iron Man 
testing: a) TRO-I anneal run means, and b) RID anneal run means. See 
Table 2 for corresponding run (0 numbers. 



The emitter aruieal cydc time was not mcxiified because the pixil 
was a 50^ improvement in total uniformity. To impmve thermal 
uni/onruty, the dynamic optimization stratcg)' was invoked during * 
the oitica) anneal step. The oxide growth during this step was com- 
puted using measured wafer temperatures in a two-stream vali- 
dated oxide growth mode) (Rg. 4). The improvement over base- 
line riD is 7x with Self-Mean calculations. The Total Thermal 
Uniformity (TTU) improvement was 18x. 

Figures 5, 6, and 7 compare measured wafer temperatures 
between MBTC and PID control for each process. In every case 
MBTC yielded much tighter control of wafer temperatures rela- 
tive to the desired setpoint. The deviation from setjx)int in the 
MBTC emitter anneal process is an expected result of tfie Dynamic 
Optimization strateg)'. Once the MBTC version of each process 
was optimized, PIXIs were performed to evaluate performance 
(see Table 1). All of the JDA goals were met or exceeded in im- 
plementing MBTC 

Measured oxide thickness uniformity and throughput were 
equal or better with MBTC. In particular, the throughput im- 
provements for the two oxidation processes correspond to ap- 
proximately 157f lower cost of ownership (calculated as cost /good 
wafer/year using the Sematech COO model). 

MBTC reliability testing (Iron Man) 

Once the capabilities of MBTC were proven, an Iron Man was per- 
formed to challenge the reliability' and robustness of the system. 
The Iron Man is Motorola's standard procedure for testing un* 
proven or recently developed, semiconductor capital equip- 
ment In Phase 1 of die Iron Man testing, various process and equxp- 
ment changes (Table 2) were made to the system to strain the 





Table 2, CquipmMt and procMS cbangas 
mada during tha PnMron Man taaling 


Rvi 


ilD# 


Equipmant changes 


1 




Profile t/c removal and replacement 


2 




Quartz process tube change 


3 




Spike t/c removal and replacement 


A 




Quartz pedestal change 


5 


iilD# 


Power acquisition module inversion 
Proceee changes 


6 




Thermocouple amplifier PCB change 


7 


< 


Gas flow change (4100% of nominal) 


8 




Gas flow change (60% of nominal) 


9 




Increased wafer load by 24 


10 




Increased water load by 48 


11 




Increased wafer load by 74 



MBTC hardware and software and detenmine failure modes. In 
Phase II, the VTR was run for an extended period of time without 
any changes to test the reliability of the entire system and simu- 
late the manufacturing environment. 

After each change was introduced, the furnace was tested 
with both MBTC and PID processes using the t/c instrumented 
wafer load. Wafer temperatures were analyzed to evaluate the 
effect of the equipment or process change on each temperature 
control method. Throughout the Phase I Iron Man testing, MBTC 

conmuedonpage I3i 
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- Furnace M9TC conthued from page 128 
performed reliabJy and proved to be less sensitive to equipment 
changes than PID. PID control degraded suffidently after the quartz 
pedestal change to cause an abort for profile t /c tenxperaturc limits 
"out of tolerance" which was set at ±5.0°C. The system had to be 
reprofiled after this change to bring PID control back into specifi- 
cation. In^iddition to the abort PID control exhibited significant 
run-to-run variation compared to MBTC. Rgure 8 compares 
predicted oxide means calculated from the wafer temi>erature 
measurements made during the Phase I Iron Man testing. The pre- 
dicted oxide means include nteasurements from all 5 t/c wafers. 

Ehmng the Phase D Iron Man, the primary objective was to eval- 
uate the reliability of the new MBTC hardware and software. 
The final system configuration was tested for 20 days under typical 
production line conditions. The laboratory furnace was run 




•iement and process chamber. 



24 hours/day, six days/ week, utilizing all three processes and both 
MBTC and PID contn>I. During idle periods, the MBTC hardware 
and software was still operational, maintaining furnace standby 
temperatures. Over the course of the Fliase n Iron Man, 1 22 jmocess 
runs were completed. The MBTC hardware and software expe- 
rienced no failures and did not cause any system douTitime. 

Conclusion 

MBTC provides precision wafer temperature control in an ox- 
idation/diffusion or LPCVD furnace. By using a thermal model 
of the wafers and wafer environment, MBTC can control 
actual wafer temperatures instead of relying on indirect meas- 
urements from thermocouples positioned in the furnace element 
and process chamber. Beyond simply controlling the furnace 
element temperature, MBTC can precisely predict and control 
the product temperature. 

We have successfully implemented MBTCona VTO7000.MBTC 
has improved silicon wafer thermal uniformity' as much as 7x and 
improved oxide film thickness unifornuties by more than 50%^ 
while achieving a 15% increase in throughput and corresponding 
reduction in cost of ownership. The reliability and production 
worthiness of the new system was demonstrated by the suc- 
cessful completion of an Iron Man test. ■ 
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